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SUMMARY 


Problem 

Devise  a method  of  calibrating  acoustic  hydrophone  line  array  segments 
that  are  about  100  meters  long. 


Results 


A calibration  method  applicable  for  long  line  array  segments  is  investi- 
gated. 

A theoretical  discussion  is  given  on  the  corrections  required  to  modify 
measured  receiver  sensitivity  curves  obtained  by  use  of  a simple  calibration 
device.  These  corrections  are  for  the  effects  of  phase  difference,  spherical 
spreading,  and  image  interference.  Curves  are  presented  from  which  the 
calibration  corrections  can  be  obtained  for  array  segments  of  3 to  24  layers 
and  for  calibration  units  h ung  radii  to  3 meters  and  covering  a frequency 
spectrum  up  to  10  kHz. 

A calibration  unit  constructed  at  NUC  in  1972  is  described,  and  test 
results  using  this  unit  are  presented.  Errors  associated  with  this  unit  due  to 
phase  differences  are  less  than  -0. 1 dB  below  1 kHz,  spherical  spreading  losses 
are  no  more  than  -0.2  dB,  and  image  interference  effects  range  from  0 dB  at 
10  Hz  to  0.1  dB  at  1 kHz,  indicating  that  no  correction  to  the  measured 
receiver  sensitivity  is  required  for  this  calibration  unit. 


Recommendation 

This  type  of  calibration  unit  should  be  used  for  the  calibration  of  all 
future  Navy  line  array  segments. 
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INTRODUCTION 

During  the  NUC  Forward  Area  Sonar  Research  cruise  in  1969  (FASOR 
III),  a short  line  array  having  an  active  hydrophone  section  of  1.7  meters  was 
used  for  acoustic  propagation  studies.  This  short  array  was  calibrated  at  Trans- 
dec  (the  NUC  test  calibration  facility)  by  suspending  the  active  section  below  a 
submerged  beam  within  the  test  facility  and  using  the  substitution  method  of 
calibration. 

This  technique  worked  for  a short  array,  but  could  not  be  used  for  an 
array  segment  over  100  meters  long.  In  fact,  a long  array  segment  had  never 
been  calibrated. 

In  December  1971,  the  Naval  Undersea  Center  (NUC)  conducted  a 
series  of  successful  acoustic  noise  tests  in  Exuma  Sound,  the  Bahamas,  using  a 
leased  long  linear  hydrophone  array.  The  results  of  these  tests  were  instrumen- 
tal in  establishing  the  Large  Aperture  Marine  Basic  Data  Array  (LAMBDA) 
project  at  NUC,  and  this  introduced  NUC  into  the  field  of  long  line  arrays. 

The  task  then  was  to  devise  a method  for  calibrating  long  line  array 
segments. 

In  January  1972,  NUC  designed  and  constructed  a calibration  unit  for 
hydrophone  line  arrays.  It  is  the  purpose  of  this  report  to  describe  and  evaluate 
this  unit  for  present  and  future  users:  for  present  users,  to  apprise  them  of 
corrections  they  should  make  to  their  calibrated  receiving  sensitivity  curves, 
and  for  future  users,  to  inform  them  of  the  unit’s  existence  and,  more  impor- 
tant, its  capabilities. 

The  report  is  divided  basically  into  two  parts.  In  the  first  part  there  will 
be  a general  discussion  of  the  method  of  calibrating  hydrophone  line  arrays.  In 
the  second  part,  NUC’s  calibration  unit  will  be  described  and  evaluated  on  the 
basis  of  its  test  results. 
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LINE  ARRAY  CALIBRATION  METHOD 
Calibration  Requirements 

The  history  of  hydrophone  calibration  methods  is  a history  of  facilities 
capable  of  performing  the  various  types  of  calibrations  that  have  been  applied 
to  hydrophones  to  produce  a value  called  the  receiving  response,  or  sensitivity, 
of  the  hydrophone. 

The  various  methods  of  calibration,  whether  absolute,  reciprocal,  sub- 
stitution, or  the  like,  furnish  a value  of  sensitivity,  representing  the  characteris- 
tics of  the  hydrophone  being  tested.  However,  these  methods,  accepted  by  the 
scientific  community  for  omnidirectional  hydrophones,  fail  when  applied  to  the 
receiving  response  of  a directional  line  array. 

A line  array  is  defined  as  a group  of  hydrophones  in  a one-dimensional 
line.  Other  lactors,  such  as  number  of  hydrophones  within  a group,  spacing 
between  hydrophones  or  between  groups,  and  operational  frequency  range, 
characterize  the  line  array  under  consideration. 

The  line  array  can  be  either  linear  or  tapered.  Linear  arrays  have  repeat- 
able  spacings  between  hydrophone  groups  or  elements,  while  tapered  arrays 
have  a spatially  different  separation  between  elements  in  the  hydrophone  group. 
In  a tapered  array  a chosen  point  within  the  group  represents  the  phase  center 
of  the  group;  that  is,  an  acoustic  wave  traveling  from  either  end  of  the  group 
will  arrive  at  the  acoustic  phase  center  in  phase. 

The  d fficulty  in  calibrating  a line  array  is  primarily  with  arranging  the 
length  of  the  array  and  the  separation  between  the  array  and  the  transmitting 
source  such  that  the  acoustic  signal  from  the  transmitter  arrives  in  phase  at  all 
hydrophone  group  centers  simultaneously. 

At  high  frequencies  f 10  kHz  or  greater),  this  problem  becomes  one  of 
having  a large  enough  tank  to  perform  the  experiment,  as  seen  in  Table  1 . 

Table  1 . Far-field  distances  as  a function  of  frequency. 


Frequency,  Hz  Range,  m 


1 

3.75  X 

I04 

10 

3.75  X 

I03 

100 

3.75  X 

102 

IK 

3.75  X 

101 

10K 

3.75 

IOOK 

3.75  X 

10' 

j 
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At  low  frequencies,  the  problem  becomes  one  of  performing  the  experi- 
ment over  many  miles  of  lake  or  ocean  where  the  receiving  array  and  transmit- 
ting transducer  must  be  at  the  proper  depth  in  isovelocity  water  to  receive  an 
acoustic  signal  that  has  not  been  interfered  with  by  multipath  arrivals. 

Before  the  introduction  of  the  NUC  calibration  unit,  the  array  had  to 
be  placed  in  the  far  field  of  the  acoustic  source.  This  range  is  given  as 

R>L:/\  (I) 

where  L is  the  length  of  the  array  and  X is  the  wavelength  of  the  transmitted 
signal. 

This  equation  can  be  expressed  as 


•»  f2 
1 o 

for  a linear  array  having  n element  groups  with  spacings  of  X 0I2 , from  which 
we  have  used  f q = c/\0  as  the  half  wavelength  frequency,  c as  the  velocity  of 
sound  in  water,  and  / as  a variable  frequency. 

This  equation  was  used  to  devise  Table  1 with  n = 1 l,c  ' * 1500  m/s 
and  f=fo- 

It  is  seen  from  Table  1 that  at  low  frequencies  (f  < 100  Hz),  line  arrays 
could  not  be  calibrated  with  existing  facilities. 


Design  Considerations 

A line  array  coiled  in  a circle  would  have  all  hydrophones  at  an  equal 
distance  from  a transmitter  located  at  the  center  of  the  coil.  This  configuration 
allows  all  signals  to  arrive  in  phase  at  all  hydrophones. 

Line  array  segments,  which  make  up  the  total  line  array,  are  usually 
around  50  to  100  meters  long.  This  large  size  necessitates  that  the  array  be 
coiled  in  a helix  and  the  problem  then  is  to  construct  this  helix  so  that  the 
measured  hydrophone  sensitivity  will  not  be  adversely  affected  by  phase  differ- 
ences, spherical  spreading,  and  image  interference. 

Phase  Differences.  For  a line  array  coiled  in  a helix,  partially  shown  in 
Fig.  1,  the  output  from  the  array  is  given  by: 

L M 
e=/  m-i 


(3) 


where  we  have  L layers  each  having  M hydrophones.  I'W|V  is  the  response  ol 
hydrophone  m in  layer  8,  and  6„lk  is  the  phase  delay  between  these  hydrophones 
For  uniform  hydrophone  response.  V,  we  have: 


where 

d is  the  diameter  of  the  line  array  hose  body 
R is  the  radial  distance  from  the  center  of  the  helix  to  the  array 

The  error  in  sensitivity  due  to  phase  differences  between  hydrophones 
in  the  different  layers  within  the  helix  can  be  determined  from: 

Phase  error  = 20  log  { VQ  /\1  V ( I + 2K) } (b ) 


(I  +2K) 


8 = J,  5,  7, ...  L (7) 


for  an  odd  number  of  layers. 


Phase  error  * 2d  log  | \ u JMl  ( I * A/}  1 


l*>J4  ( / * cos  4 ) N sin  4^  ) 

\ Lift  / \ft  /J  // 


n * k> 


4 4.  /».  X. 


lor  an  even  number  o!  layers. 


where 


- 2 


K=.1  or  4 


V - running  index  lo  account  lor  the  even  or  odd  / hydrophones 


Computer  programs  were  written  to  evaluate  these  and  the  follow  mg 
equations,  and  are  listed  in  the  Appendix. 

IX*sign  considerations  have  been  given  lo  include  frequencies  and  num- 
ber of  layers  far  beyond  the  practical  limit  to  show  the  onset  ol  phase  inter- 
ference. 

Figures  2 through  I are  presented  to  show  the  phase-  errors  at  I . 5.  and 
10  kHz  versus  helix  radii  from  ().<»  to  3 meters  and  lor  helixes  trom  3 to  I ~ and  3 
to  24  layers  and  an  array  hose  diameter  of  l>  cm.  At  high  frequencies  ol  5 or 
10  kHz.  an  interference  pattern  is  established  lor  a large  number  ol  layers,  and 
the  phase  error  oscillates  about  the  -10  dB  level.  Lay  ers  greater  than  I 7 were 
therefore  omitted  to  avoid  confusion  of  the  sort  seen  in  I ig.  4 

Phase  errors  were  found  to  decrease  with  increasing  radii. 


Spherical  Spreading.  Spherical  spreading  losses  due  to  the  difference  in 
separation  between  source  transducer  and  reference  hydrophone  and  source  and 
array  segment  layer  were  derived  similarly  and  are 


Distance  error  = Jl)  log 


"-"Ik 


'( i * :ki\ 


V = 5. 


for  an  odd  number  of  layers. 


I 


Distance  error  = 20  log 


(I  + K) 


»**-U j 

<i=-4  7 


If  = •/.  ft.  .v z. 


(12) 


for  an  even  number  of  layers, 
where 


(13) 

(14) 


Figure  5 illustrates  the  error  in  spherical  spreading  for  helix  radii  from 
0.6  to  3 meters  and  for  a helix  of  3 to  24  layers.  Spherical  spreading  er.ors 
were  found  to  decrease  with  increasing  mdii  and  to  increase  with  incr<  asing 
number  of  layers. 


Image  Interference.  Image  interference,  or  Lloyd  mirror  effect,  must  be 
taken  into  consideration  whenever  an  acoustic  source  and  receiver  are  near  the 
surface  or  the  bottom.  This  is  due  to  constructive  and  destructive  interference 
between  the  direct  signal,  the  surface  retlected  signal,  and  the  bottom  reflected 
signal. 

If  we  extend  the  theory  as  given  by  Urick*  to  include  bottom  reflected 
sound  and  spherical  spreading,  we  can  generate  another  family  of  curves  to 
show  the  losses  due  to  image  interference. 

Let 


A . , 

= — sin  f c ot) 

Rs 

(15) 

A . . 

P,  =-P/  £-sin  u(t  + T,j) 

(16) 

A . , 

p,  = P-,  — sin  c o(t  + 7,  J 

J ‘ A-. 

(17) 

Urick, R.J.  Principles  of  underwater  sound  for  engineers,  pp.  110-13.  McGraw-Hill  Book  Co.,  1967. 


where  P3.P1,  and  p2  represent  ihe  pressures  at  a field  point  for  the  direct  signal 
and  the  surface  and  bottom  reflected  signals,  respectively  (Fig.  6). 

A is  the  source  strength  at  the  source  point 

p,  and  p2  are  the  surface  and  bottom  reflection  coefficients 

tij=(R,  -R3)/c  (18) 

t23  ~ (^2  ~ ^3^^  (19) 

t = time  differences  between  the  different  arrivals 

The  combined  pressure  at  a field  point  on  one  layer  is  given  by 

P=P,+P2+P3  (20) 

and  the  intensity  is  given  by  the  time  average  of  the  squared  pressure  over  a half 
cycle  as 


(P,  + P2  +P3>2 


» Tl  ' 


+p2  +p3)2  dt 


where  p = the  density  of  water. 

After  performing  the  integration  and  introducing  the  average  intensity 
at  the  field  point  due  to  the  source  alone  as 


P3  A2 
° 2 Pc  2pcR32 


we  have 


/RjV  /R,\2 

g +2Ri 


H2  COs(<jJT2j)  Hi  COS( UTjj) 


R , 


2 /?/ 

- — — 1 jcosf'tJT/i>  cos (ut,3)  + s\n(uiTl3)  s\n(uT23)\ 

KlK2 


(25) 


Equation  25  thus  shows  the  modification  of  the  acoustic  field  due  to  image 
interference  of  the  surface  and  bottom  at  a field  point  on  one  layer. 

The  total  image  error  for  all  layers  of  the  array  segment  is  then 


Image  error  = 10  log 


= 10  log 


1 = 2,  4.6,...  L 


1 = 3,  5,  7 L 


(26) 


(27) 


Figures  7 through  1 1 show  the  changes  in  image  error  versus  radii  for  the 
various  even  and  odd  layers  at  1,5,  and  10  kHz.  These  values  were  obtained  by 
considering  all  reflection  coefficients  as  n = 1 , and  the  calibration  unit  midway 
between  surface  and  bottom.  In  all  cases  the  trend  is  to  increase  image  errors 
with  increasing  radii  and  increasing  error  for  increasing  number  of  layers  at 
1 kHz  up  to  1 2 or  1 3 layers,  at  which  time  the  error  decreases.  At  5 kHz  the 
error  curves  are  fairly  well  grouped  together,  and  at  10  kHz  the  curves  diverge 
into  two  families  of  curves. 
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Calibration  Corrections.  In  each  instance  of  phase  differences,  radial 
spreading,  and  image  interference,  the  value  obtained  for  the  measured  receiving 
sensitivity  of  the  array  segment  will  have  been  reduced  for  all  negative  error 
terms  and  increased  for  those  that  were  positive.  To  remedy  this  the  corrected 
array  receiver  sensitivity  for  small  errors  is  given  algebraically  by: 

Corrected  array  sensitivity  = Measured  receiver  sensitivity 

- Phase  error  - Distance  error 

- Image  error  (28) 
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Measured  receiver  sensitivity  of  a hydrophone  is  obtained  by  the  substi- 
tution method  of  calibration  and  is  usually  modified  by  image  interference 
effects  between  the  source  transducer,  the  calibrated  reference  hydrophone, 
the  water  surface,  and  the  bottom.  For  frequencies  below  1 kHz  and  for  near- 
normal incidence  from  surface  and  bottom,  the  assumption  /a  = / is  quite  valid. 
Under  these  circumstances,  if  the  source  and  reference  hydrophone  are  midway 
between  surface  and  bottom,  Eq.  25  shows  that  ///0  = 1 and  this  extra  modifi- 
cation does  not  need  to  be  made. 


figure  I . layer  counting,  one  side  of  helix. 
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Helix  Radius 

Distance  error  versus  helix 


NUC  CALIBRATION  UNIT 


Unit  Description 

The  NUC  hydrophone  line  array  calibration  unit  (U.  S.  Patent  No. 
3,859,620,  issued  5 February  1975)  was  constructed  at  NUC  in  January  1972 
and  is  shown  in  Fig.  1 2. 

The  sketch  in  Fig.  13  shows  the  top  and  base  (2.7-meter-diameter 
discs),  which  were  made  from  double-thickness  1 .9-cm  plywood  coated  with 
three  coats  of  fiberglass  resin.  The  top  disc  has  a 0.9-meter-diameter  hole  in  the 
center,  and  the  disc  has  been  drilled  for  four  1 .0-cm  eyebolts  which  attach  to  a 
lifting  harness.  The  top  disc  is  also  drilled  for  eight  0.6-cm  eyebolts,  from  which 
the  5.5-meter,  1 ,3-cm  line  holds  the  calibration  housing.  This  housing  is  made 
of  heavy  wire  mesh  which  is  stiff  enough  to  retain  a cylindrical  shape  during 
calibration  and  is  opaque  to  acoustic  transmissions  within  the  useful  frequency 
range. 

The  unit  is  assembled  by  fitting  the  wire  mesh  over  the  base,  which  is 
supported  from  the  ground  by  wooden  beams.  The  wire  mesh  is  then  attached 
temporarily  to  the  base  with  large,  easily  removable  staples.  The  array  is  then 
coiled  within  the  wire  mesh  frame  and  secured  with  line,  as  shown  in  Fig.  14. 

The  top  support  is  then  installed  within  the  coiled  array,  preventing  the  top 
from  pressing  on  the  array  during  shipment  to  the  calibration  facility.  The  5.5- 
meter  support  lines  are  attached  to  the  wire  mesh  housing,  the  top  disc  is  low- 
ered within  the  wire  frame,  and  staples  are  applied  to  hold  the  top  in  place 
during  shipment  (Fig.  15). 

Upon  the  unit’s  arrival  at  the  calibration  facility,  all  staples  are  removed 
and  a crane  lifts  out  the  top  disc.  The  support  cables  extend  until  the  calibra- 
tion unit  is  lifted  free  of  the  truck  bed,  as  shown  in  Fig.  1 2.  The  unit  is  then 
placed  within  the  calibration  tank  and  is  free-floated  into  position  for  calibration. 

Test  Results 

In  this  section  of  the  report  the  line  array  calibration  unit  is  evaluated 
for  errors  and  for  future  users,  and  the  results  of  the  initial  test  of  a standard 
geophysical  array  manufactured  by  the  Seismic  Engineering  Company  are  pre- 
sented. This  array  occupied  8 layers  within  the  calibration  unit. 


Phase  Differences  for  NUC  Unit.  Figures  16  through  19  exhibit  the  phase 
errors  over  the  frequency  range  of  1 00  Hz  to  1 0 kHz  for  3 to  24  layers  of  a 
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line  array  segment  within  this  helical  wire  frame.  Phase  errors  were  found  to  be 
negligible  below  100  Hz  for  all  layers,  to  increase  with  number  of  layers,  and  to 
increase  with  frequency.  For  the  geophysical  array  tested,  errors  due  to  phase 
differences  were  less  than  -0.1  dB  below  1 kHz. 


Spherical  Spreading  for  NUC  Unit.  Spherical  spreading  losses  were  seen 
to  be  0.2  dB  from  Fig.  5,  for  the  NUC  calibration  unit  and  geophysical  array. 


Image  Interference  for  NUC  Unit.  Figures  20  through  23  show  the  loss  in 
intensity  ever  the  frequency  range  of  1 00  Hz  to  1 0 kHz  for  3 to  24  layers  of  a line 
array  segment  within  the  NUC  calibration  unit.  For  this  unit,  image  errors  were 
found  to  increase  with  number  of  layers  and  with  increasing  frequency  up  to 
1 kHz.  Above  1 kHz  the  image  error  oscillates  about  0.15  dB  until  the  curves 
break  up  into  two  families  of  curves  for  the  odd  and  even  layers  near  9 kHz. 

For  the  geophysical  array  used,  image  errors  were  approximately  0 dB  up  to 
100  Hz  and  increased  to  0.1  dB  at  1 kHz. 


Calibration  Corrections  for  NUC  Unit.  Figure  24*  shows  the  results  of 
an  18  February  1972  measurement  of  the  receiver  sensitivity  for  the  standard 
Seismic  geophysical  array  from  40  Hz  to  1 kHz.  Corrections  to  this  curve  were 
found  to  be  negligible,  being  -0.1  dB  at  40  Hz  and  -0.2  dB  at  1 kHz. 


F rom  Naval  Undersea  Center  report  6480,  “Transducer  calibration  of  Seismic  streamers,"  July  1972. 
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Figure  16.  Phase  error  versus  frequency,  NIT  unit,  odd  layers,  100  Hr  - 1 kHz. 


Figure  17.  Phase  error  versus  frequency,  NIT  unil.cven  layers,  100  Hz  - 1 kHz. 
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Figure  24.  Measured  Seismic  array  sensitivity. 


SUMMARY  OF  RESULTS 

V 

1.  Phase  errors  were  found  to  be  negligible  below  100  Hz  for  all  layers, 
to  increase  with  number  of  layers,  to  increase  with  frequency,  and  to  decrease 
with  increasing  radii. 

2.  Spherical  spreading  errors  were  found  to  decrease  with  increasing 

radii. 

3.  Image  errors  were  found  to  increase  with  number  of  layers  and  with 
increasing  frequency  up  to  1 kHz.  Above  1 kHz  the  image  error  oscillates  about 
0.15  dB  until  the  curves  break  up  into  two  families  of  curves  for  the  odd  and 
even  layers  near  9 kHz.  Image  errors  also  increase  with  increasing  radii  and 
increasing  number  of  layers  at  1 kHz  up  to  1 2 or  1 3 layers,  and  then  decrease. 

At  5 kHz  the  image  error  versus  radii  curves  are  fairly  well  grouped  together, 
and  at  1 0 kHz  the  curves  diverge  into  two  families  of  curves. 
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CONCLUSIONS 


A calibration  unit  constructed  at  NUC  in  early  1972  has  been  used  to 
successfully  calibrate  a large  number  and  variety  of  line  array  segments. 

It  has  been  shown  that  errors  associated  with  this  unit  due  to  phase 
differences  between  the  layers  of  the  helix  are  less  than  -0.1  dB  below  1 kHz, 
that  spherical  spreading  losses  are  no  more  than  -0.2  dB,  and  that  image  inter- 
ference effects  range  from  0 dB  at  10  Hz  to  0.1  dB  at  1 kHz,  indicating  that  no 
correction  to  the  measured  receiver  sensitivity  is  required. 


RECOMMENDATIONS 

Considering  the  successful  calibration  measurements  that  have  been 
made  with  this  unit,  it  is  recommended  that  continued  use  be  made  of  this  unit, 
or  a newer  version  of  it,  for  all  future  NUC  line  array  segments  and  that  similar 
units  be  used  at  all  Navy  calibration  facilities. 


Appendix 

COMPUTER  PROGRAMS  FOR  HP  9830 
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1 REM  SPHERICAL  SPREADING  ERRORS 
10  DIM  Z [25, 10] 

20  D=3 . 5 

30  FOR  R=2  TO  10 
40  YO=Y8=0 
50  GOSUB  320 
60  NEXT  R 

70  SCALE  0,11, -2. 7,0.1 
80  XAXIS  -2.6,1,2,10 
90  YAXIS  2, 0.2, -2. 6,0 


100 

LABEL  (“, 2 

, 1.7,0,8/11) 

no 

FOR  X=2  TO 

10 

120 

PLOT  X-0.3 

rH 

CM 

l 

130 

CPLOT  2,-0 

.3 

140 

LABEL  C 1 50 

)X 

150 

FORMAT  F3 . 

0 

160 

NEXT  X 

170 

LABEL  ( ::,  2 

,1.7,0,8/11) 

180 

FOR  Y=- 2 . 6 

TO  0 STEP  0.2 

190 

PLOT  0 , Y,  1 

200 

CPLOT  -2,0 

. 1 

210 

LABEL  (220) Y 

220 

FORMAT  F 1 1 

. 1 

230 

NEXT  Y 

240 

FOR  L=3  TO 

24 

250 

PEN 

260 

FOR  R=2  TO 

10 

270 

PLOT  R, Z [ L 

,R] 

280 

NEXT  R 

290 

NEXT  L 

300 

LETTER 

310 

STOP 

320 

M=3 

330 

FOR  G= 1 TO 

2 

340 

K=  1 

350 

FOR  L=M  TO 

(M+2 1 ) STEP  2 

360 

A0=(L-l)"D/24 

370 

PO  = SQR(A0!\ 

A0  + R::R) 

380 

IF  (L/2-K) 

<0.7  THEN  450 

390 

Y0=Y0+1/P0 

400 

R0  = SQR((D!! 

D/576)+R"R) 

410 

A=(1+R05!Y0 

)/ ( 1+K) 

420 

IF  A> 1 THEN  500 

430 

Z[L,R]=20” 

LGT(A) 

440 

GOTO  510 
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450  Y8=Y8+1/P0 

460  B = (1  + 2::R”-Y8)/(1  + 2”K.) 

470  IF  B> 1 THEN  500 
480  Z[L,R]=20:;LGT(B) 

490  GOTO  510 
500  Z [ L, R] =0 

510  WRITE  (15,520)L,R,Z[L,R] 

520  FORMAT  2F10.1,F16.1 

530  K=K+ 1 

540  NEXT  L 

550  M=M+1 

560  NEXT  G 

570  RETURN 

580  END 


1 REM  PHASE  INTERFERENCE  ERRORS 
10  DIM  V[25] , P[25, 10] 

20  C=5000 
30  D=3  . 5 
40  R=4 . 5 
50  F= 1 0 0 

60  FOR  N=1  TO  10 

70  S0=C0=0 

80  W = 2::PI::F/C 

90  GOSUB  370 

100  F=F+100 

110  WRITE  (15,520) 

120  NEXT  N 

130  SCALE  0,1100,-1.1,0.1 
140  XAXIS  -1,100,100,1000 
150  YAXIS  100,0.1,-1,0 
160  LABEL  (”,2,1.7,0,8/11) 

170  FOR  X=100  TO  1000  STEP  100 

180  PLOT  X-100, -1.1,1 

190  CPLOT  2,-2 

200  LABEL  (280)X 

210  NEXT  X 

220  LABEL  (:!,  2 , 1 . 7 , 0 , 8 / 1 1 ) 

230  FOR  Y = - 1 TO  0 STEP  0 . 1 
240  PLOT  0 , Y , 1 
250  CPLOT  -2,0.1 

260  LABEL  (261)Y 

261  FORMAT  F4.1 
270  NEXT  Y 


280 

FORMAT  F8.0 

290 

L=  2 3 

300 

F=  100 

310 

PEN 

320 

FOR  N= 1 TO  10 

330 

PLOT  F,P[L,N] 

340 

F=F+100 

350 

NEXT  N 

351 

LETTER 

360 

STOP 

370 

K=  1 

380 

FOR  L=3  TO  23  STEP  2 

390 

A0  = (L-l)::D/24 

400 

P0=SQR(A0"A0+R::R) 

410 

X0=W5:(P0-R) 

420 

CO=CO+COSCXO) 

430 

S0=S0+SI ncxo ) 

440 

T0=SQR(C0"C0+S0::S0) 

450 

V[L]=ABS((l  + 2::T0)/Cl  + 2 

::K)) 

460 

IF  V[L]<0.0001  THEN  480 

470 

P[L,N]=20::LGT(V[L]  ) 

480 

GOTO  500 

490 

P[L,N]=-80 

500 

WRITE  Cl5,510)L,F,V[L] 

,P[L,N] 

510 

FORMAT  2F10 . 1,7X,E9 .2, 

F 1 3 . 1 

520 

FORMAT  / 

530 

K = K+ 1 

540 

NEXT  L 

550 

RETURN 

560 

END 

10 

REM  IMAGE  INTERFERENCE 

ERRORS 

20 

DIM  V [ 2 5 ] , P [ 2 5 , 1 0 ] 

30 

C = 5000 

40 

D = 3 . 5 

50 

R = 4 . 5 

60 

R9=R”R 

70 

R0=SQRCR9+1296) 

80 

10  = 1 

90 

F=  1 0 0 

100 

FOR  N= 1 TO  10 

1 10 

S0=C0=I 1=0 

120 

W=2”P  I :tF/C 

130 

GOSUB  440 

1 


140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 


F=F+100 

WRITE  (15,710) 

NEXT  N 

SCALE  0,1100,-0.05,0.32 
XAXIS  0,100,100,1000 
YAXIS  100,0.05,0,0.3 
LABEL  (5:, 2, 1.7, 0,8/11) 

FOR  X=100  TO  1000  STEP  100 

PLOT  X-40,0,1 

CPLOT  0,-1 

LABEL  (330)X 

NEXT  X 

LABEL  (”,2,1.7,0,8/11) 

FOR  Y=0  TO  0.3  STEP  0.05 

PLOT  0 , Y, 1 

CPLOT  -4,0.1 

LABEL  ( 3 1 0 ) Y 

FORMAT  F6.2 

NEXT  Y 

FORMAT  F5.0 

FOR  L = 4 TO  24  STEP  2 

F=  1 0 0 

PEN 

FOR  N=1  TO  10 
PLOT  F, P [L, N] 

F=F+100 
NEXT  N 
NEXT  L 
LETTER 
STOP 
S = 2 

FOR  Q=1  TO  2 
11  = 0 

FOR  L=S  TO  (S+23)  STEP  2 

Ml=(L-l)/2 

M=-M1"D/12 

R7  = R9+M::M 

A=SQR(R7 ) 

FOR  G= 1 TO  2 

Rl  = SQR(R9+(36-M)::(36-M)) 

R2  = SQR(R9  + (36+M)::(36+M)) 

T 1 =( R 1 -A) 

T2=(R2-A) 

U = ((2::R7)/(R1::R2))::(C0S(W::T1)::C0S(W”T2) 
+ SIN(W!!T1)5!SIN(W::T2)) 


29 


58  0 I1  = I1  + 1 + R7”(1/(R1::R1)  + 1/  (R2::R2)  ) + 2::A:c(COS( W::T2 ) 
/R2-C0S(W::T1)/R1)-U 
590  M=-M 
600  NEXT  G 
610  IF  Q= 1 THEN  640 
620  V[L] =ABSCC I 1+1 0 ) / L) 

630  GOTO  650 
640  V [ L] =ABS ( 1 1 / L) 

650  IF  V[L]<0.0001  THEN  680 
660  P[L,N]  = 10::LGT(V[L]) 

670  GOTO  690 
680  P[L,N]=-80 

690  WRITE  (15,700)L,F,V[L],P[L,N],R1,R2,A/U/I1 
700  FORMAT  2F8.1,3X,E9.2,F10.1,3F8.1,2F8.3 
710  FORMAT  / 

720  NEXT  L 
730  S=S+1 
740  NEXT  Q 
750  RETURN 
760  END 
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